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ABSTRACT 

Chandra ACIS observations of PKS 0521— 365 find that the X-ray emission of this 
BL Lac object consists of emission from an unresolved core, a diffuse halo, and a 
2 arcsec jet feature coincident with the inner radio/optical jet. A comparison with a 
new ATCA 8.6-GHz map also finds X-ray emission from the bright hotspot south-east 
of the nucleus. The jet spectrum, from radio to X-ray, is probably synchrotron emission 
from an electron population with a broken power-law energy distribution, and resem- 
bles the spectra seen from the jets of low-power (FR I) radio galaxies. The hotspot 
X-ray flux is consistent with the expectations of synchrotron self-Compton emission 
from a plasma close to equipartition, as seen in studies of high-power (FR II) radio 
galaxies. While the angular structure of the halo is similar to that found by an analysis 
of the ROSAT HRI image, its brightness is seen to be lower with Chandra, and the 
halo is best interpreted as thermal emission from an atmosphere of similar luminosity 
to the halos around FR I radio galaxies. The X-ray properties of PKS 0521— 365 are 
consistent with it being a foreshortened, beamed, radio galaxy. 

Key words: galaxies: active - BL Lacertae objects: individual: PKS 0521— 365 - 
X-rays: galaxies - radio continuum: galaxies 



It has been suggested (e.g., Antonucci & Ulvestad 1985; 
Ulrich 1989; Urry & Padovani 1995) that BL Lac objects 
are "unified" with FR I radio galaxies (FRIRGs), so that 
the only difference between BL Lac objects and FRIRGs 
is their orientation relative to the observer's line of sight. 
Thus BL Lac objects are supposed to be viewed from a 
direction close to the nuclear jet, whose Doppler-boosted 
emission dominates the appearance of the object at radio to 
7-ray energies. FRIRGs are supposed to be viewed from a 
direction outside the nuclear jet, so that the nucleus is faint 
relative to the outer radio structure, the optical galaxy, and 
the unbeamed X-ray emission. 

This idea has been tested through comparisons of 
BL Lac and FRIRG number densities, extended (and hence 
isotropic) radio luminosities and morphologies, and the opti- 
cal luminosities and morphologies of the underlying galaxies 
(Morris et al. 1991; Perlman & Stocke 1993; Pesce, Falomo 
& Treves 1996). Tests based on X-ray properties, which are 
intrinsically among the most interesting since strong X-ray 
emission is a defining feature of BL Lac objects, have fitted 
the distribution of X-ray luminosity to a theoretical function 
based on some beaming factor and beam opening angle, us- 
ing either the total X-ray luminosity (e.g., Padovani & Urry 
1990) or the core X-ray luminosity, after any extended X-ray 
emission has been removed (Canosa 2000; Hardcastle et al. 
2001). 



Another possible X-ray test of the identification of 
BL Lac objects as beamed FRIRGs can be made based on 
the association of FRIRGs with strong X-ray emission from 
galaxy-to-group scale gas distributions (Worrall & Birkin- 
shaw 1994, 2000). Since this emission is thermal and ex- 
tended, it must radiate isotropically, and therefore should 
appear with similar properties around BL Lac objects and 
FRIRGs. A clear difference between the extended X-ray 
emissions near BL Lacs and FRIRGs would constitute evi- 
dence against a strong unification hypothesis. 

This test is difficult because the extended X-ray emis- 
sion has low contrast relative to the BL Lac's unresolved 
active galactic nucleus (AGN) . Wings on the point response 
function (PRF) of the X-ray optics and detector used to ob- 
serve a BL Lac are therefore likely to exceed the brightness of 
the extended emission. Nevertheless, Hardcastle, Worrall & 
Birkinshaw (1999), using the ROSAT HRI, reported a halo 
of resolved emission around the AGN in the low-redshift 
= 0.05534; Keel 1985) BL Lac PKS 0521-365. Although 
PKS 0521—365 had a low X-ray flux during this observation, 
the halo was only detectable above the wings on the HRI's 
PRF over a small range of angles. The inferred halo lumi- 
nosity, of L x (0.2 - 1.9 keV) w 8 x 10 35 W, was far higher 
than the values < 10 35 W typical of the FRIRGs in the B2 
sample of Worrall & Birkinshaw (1994). Since the core ra- 
dius of the PKS 0521— 365 X-ray halo was found to be only 
8 arcsec (12 kpc for Ho = 50 kms _1 kpc _1 , qo = 0, as as- 
sumed throughout this paper), the hot gas appeared to be 
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confined near the host galaxy, rather than being associated 
with the group of galaxies in which PKS 0521— 365 lies, as 
is more usual for FRIRGs.. 

At this luminosity the cooling time of gas in the X-ray 
halo is less than 3 x 10 s yr, so that a halo of this size and 
luminosity would be expected to be participating in a rapid 
cooling flow, with an estimated infall rate > 20 Mq year -1 . 
The ROSATHRl data could not measure the temperature of 
the gas, or say whether the extended halo had the structure 
of a cooling flow near the AGN. 

For these reasons, we obtained an imaging observation 
of PKS 0521-365 with the Chandra ACIS-S. We anticipated 
that this would provide detailed structural and spectral in- 
formation on the various components of the source, including 
the ROSAT-detected halo. 

Chandra has shown that X-ray jets are relatively com- 
mon in low-power radio galaxies (Worrall, Birkinshaw & 
Hardcastle 2001). Such jets are best interpreted as syn- 
chrotron emission from a high-energy continuation of the 
same electron spectrum responsible for the radio output. 
Thus if BL Lac objects are simply FRIRGs seen in pro- 
jection, we might expect to see the same X-radiation from 
the jets as is seen in radio galaxies, with a similar X-ray to 
radio power ratio. However, this requires that the observa- 
tion has enough angular resolution to measure a foreshort- 
ened jet feature close to the Doppler-enhanced X-ray core. 
PKS 0521— 365 is known to exhibit a radio/optical jet on 
arcsecond scales (Keel 1986; Macchetto et al. 1991) which 
connects to an inner VLBI jet on milliarcsec scales (Tingay 
et al. 1996), and so the arcsec jet might be resolvable from 
the bright core by Chandra. 

Our observation of PKS 0521— 365, made early in the 
Chandra mission, successfully detected and studied the halo, 
and found X-radiation associated with both the 1.5 arcsec- 
scale inner jet of the BL Lac and the hotspot of its radio 
structure. In this paper we describe these Chandra results, 
and their interpretation in terms of emission mechanisms 
and unified schemes for BL Lac objects. 



2 OBSERVATIONS 
2.1 X-ray observations 

The Chandra observation of PKS 0521— 365 was made on 
31 December 1999. The telescope was pointed so that 
the source was located near the aim-point of the back- 
illuminated CCD chip S3 of the Advanced CCD Imaging 
Spectrometer (ACIS). Unfortunately, at this early stage 
of the mission the standard offsets placed PKS 0521— 365 
close to the boundary between two readout nodes of the S3 
chip, so that the source moved between regions of the chip 
with somewhat different responses, and the dead region be- 
tween them, as the pointing direction of the observatory was 
dithered by about 20 arcsec. 

Since PKS 0521—365 was expected to produce a high 
X-ray count rate, in excess of 1 count s _1 , the observation 
was made with the minimum permitted ACIS frame time 
(0.4 sec). To achieve this it was necessary to observe with a 
small window on the CCD. This window was 128 pixels (each 
of 0.492 arcsec on a side) wide, and 1024 pixels long, and 
at the roll of the satellite during these observations (340°.4), 
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Figure 1. The central part of the subsampled (0.0984 arcsec 
square pixels), exposure corrected, 0.3 — 7.0 keV Chandra ACIS- 
S image of PKS 0521— 365. The linear greyscale corresponds to 
to 4 X 10 — 7 photons cm -2 s — 1 pixel -1 . The radio image super- 
imposed is our 8.6-GHz map made with the Australia Telescope 
Compact Array with 1 arcsec resolution, contoured at intervals 
of 20 mjy from 20 to 200 mjy. 



the long axis of the viewing window lies in position angle 
-70° 

Even with the reduced frame time, there is a signifi- 
cant chance that several photons will be received from the 
core of PKS 0521— 365 in a single CCD pixel within a sin- 
gle integration. This "pileup" effect leads to a distortion of 
the spectrum of the core and a reduction in the count rate 
relative to the true incoming rate. Another consequence of 
pileup is that the image of the core differs from the images of 
fainter, unresolved, sources, in the sense of appearing some- 
what extended and flat-topped. Finally, because the count 
rate is high, there is a significant probability that the CCD 
will record photons from the source while the data are being 
clocked to the readout registers. This causes faint "readout 
streaks" , which appear in position angles 20° and 200° from 
the core of the source for the —70° roll of these data. 

The data were provided to us in a variety of processing 
states: the analysis reported here was performed on data 
processed using pipeline software version R4CU5UPD11.1, 
and used software from the CIAO v2.0.2. After screening 
out periods of high background, the original exposure of 
9122 sec (live time, from an on-source time of 9870 sec) was 
reduced to 7776 sec of live time. 

An image was made from the events file by restrict- 
ing the apparent photon energies to the relatively well- 
calibrated range 0.3 — 7.0 keV, and corrections were made 
for the effective area of the instrument and the exposure. 
At the same time, the data were reprocessed to remove the 
anti-aliasing pixel randomization included in the standard 
pipeline, to produce an image with the best possible an- 
gular resolution, and were resampled onto a pixel scale one 
fifth that of the ACIS-S pixels (0.0984 arcsec). The resulting 
image is shown in Fig. oj. Clearly visible in this image is the 
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Figure 2. A 2.5-GHz radio map of PKS 0521-365 made with 
the Australia Telescope Compact Array, with 4 arcsec resolution, 
and contoured at intervals of 20 mjy from 20 to 100 mjy, then 
every 200 mjy. Bright compact components are associated with 
the core and the hotspot to its south-east. A halo of lower surface 
brightness emission surrounds these components. 



bright concentration of counts at the core of the AGN, an 
asymmetric excess of counts associated with the brightest 
feature in the radio jet to the north-west (at J2000 coor- 
dinates RA 05 h 22 m 57 ! !86, Dec -36°27'30".l), and a radial 
gradient in the count density from the AGN. An excess of 
counts is also be associated with the bright hotspot of the 
radio structure, to the south-west of the source. 

Spectral analyses were performed using the CIAO 
SHERPA package, after calculating the response matrices 
in the standard fashion. Although the matrices used for fit- 
ting were appropriate for only one of the readout nodes of 
the ACIS-S3 chip (the node on which most of the counts 
lie), tests using matrices calculated for the other readout 
node showed no significant differences from the results pre- 
sented here. Fits were performed only in the relatively well- 
calibrated energy range 0.3 — 7.0 keV, and only after the 
counts had been binned such that there were 20 or more 
counts in each energy bin. 

A radial profile centred on the core of PKS 0521— 365 
was extracted from the subsampled, 0.3 — 7.0 keV image, and 
analysed in the manner described in Worrall et al. (2001), us- 
ing an averaged point response function based on the CIAO 
PRF library and the spectra of several similar sources. This 
radial profile (shown in Fig. |^) is expected to be a poor 
representation of the true source structure near the core be- 
cause of pileup: fits to the profile therefore ignored the cen- 
tral 1.2 arcsec diameter circle. In addition, the profile was 
constructed after omitting the jet region and the readout 
streaks. The final profile contains only about a third of the 
detected counts within a radius of 23.5 arcsec. 
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Figure 3. The radial profile of PKS 0521-365 extracted from the 
subsampled image (Fig. |l]), shown as datapoints. Superimposed is 
the best-fitting composite model made up of a point source (using 
a point response function model for the Chandra ACIS-S of the 
form discussed by Worrall et al. 2001) and a beta-model with 
(3 = 0.90 and the best-fitting core radius of 9.3 arcsec. The fit in 
the inner 2 arcsec is poor, most likely because of the effects of the 
readout streak and pileup, but the outer region is well described 
by a superposed beta model and the wings of the point source. 



2.2 Radio observations 

To complement these X-ray data, we mapped PKS 0521—365 
using the Australia Telescope Compact Array (ATCA) in its 
6D configuration on 4 and 5 April 2000, giving baselines from 
77 to 5878 m. Data were obtained at 1.384, 2.496, 4.800, 
and 8.640 GHz with 128 MHz bandwidth subdivided into 
32 channels (of which 16 are independent and only 13 are 
useful). 10-second integrations were used at all bands. Ob- 
servations of PKS 0521—365 were calibrated for amplitude 
against the standard southern calibrator PKS 1934—638, 
and PKS 0537-441 was used as the local phase calibrator. 

All six antennas of the ATCA were available for most 
of the observation, although difficulties with the correlator 
caused some uv coverage to be lost, and some of the data 
were taken in bad weather, so that the phase stability was 
sometimes poor. Phase calibrations were made every 15 min- 
utes to compensate for the unstable conditions, but after 
basic calibration and data-editing using MIRIAD, a lengthy 
process of self-calibration (in MPS) was needed to correct 
the phases and improve the images. Only a small fraction of 
the data were lost because of bad amplitudes, interference, 
or phase jumps. Residual phase errors limit the dynamic 
range of the final image to about 10 4 . 

Fig. |l| shows contours of the 8.6-GHz, 1-arcsec resolu- 
tion radio map of PKS 0521—365 superimposed on an image 
of the 0.3-7.0 keV X-ray counts derived from our Chandra 
ACIS-S data. The radio map shows an unresolved compo- 
nent at the core of PKS 0521— 365, a clear jet to the north- 
west dominated by a knot about 1.5 arcsec from the core, 
and a compact hotspot to the south-east. A Gaussian fit to 
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Figure 4. The spectrum of the core of PKS 0521— 365 in en- 
ergy range 0.3 — 7.0 keV, for which the calibration is expected to 
be adequate, with the best-fitting absorbed power-law spectrum 
(upper panel). The lower panel shows the residuals from this fit. 
Note the systematic misfit near 2 keV, which arises from pileup 
of photons near 1 keV, where the count rate peaks. The fit is of 
adequate quality (\ 2 = 215 with 199 degrees of freedom), but 
biased by pilcup. 



the hotspot measures a FW HM of 0.4 arcsec. This result 
will be important in Sec. 



3.4 



Much low-surface-brightness emission does not appear 
on Fig. |l| at these contour levels, but is revealed on a lower- 
resolution map of PKS 0521-365, made with our ATCA 2.5- 
GHz data, and shown in Fig. |[ The diffuse, extended, radio 
emission associated with the source is seen to surround the 
radio core and the hotspot to the south-east. Only a slight 
extension of the radio core to the north-west betrays the 
presence of the radio jet. 

The distorted radio morphology of PKS 0521— 365 evi- 
dent in these maps is characteristic of radio-selected BL Lac 
objects, but PKS 0521—365 is somewhat unusual in showing 
the bright jet typical of FR I radio galaxies to the north-west 
combined with a hotspot, typical of FR II radio galaxies, 
to the south-east. FR II-like and FR I-like extended radio 
morphologies are about equally common in radio-selected 
BL Lac objects (Rector & Stocke 2001). It is not implausible, 
therefore, that a BL Lac should show a mixed-morphology 
radio structure, which may be related to the structures of 
radio galaxies with luminosities near the FR I/FR II break 
(Gopal Krishna & Wiita 2000). 



3 ANALYSIS 
3.1 The core 

The X-ray image of PKS 0521-365 in Fig. | is domi- 
nated by a strong core component, responsible for 1.01 ± 
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Figure 5. x 2 contours for fits to the spectrum of PKS 0521— 365 
as a function of the photon index of the fitted spectrum, T, and 
the intrinsic absorbing column at the core of PKS 0521— 365, Nn. 
Contours are drawn at 1, 2, and 3<r uncertainty offsets from the 
best fit. No excess Nn at PKS 0521— 365 is required by these fits. 



0.02 counts -1 in a 2-arcsec diameter circle. At this count 
rate the core is strongly piled up, and individual counts fre- 
quently correspond to the combined energies of two or three 
photons that arrived during a single frame. Nevertheless a 
spectrum was extracted and a spectral fit was attempted, 
using counts from a region 2.2 arcsec in diameter about the 
centre of the core and the standard choice of event grades 
(ASCA grades 0, 2, 3, 4, and 6). The effects of pileup in the 
resulting spectrum (Fig. ^|) are clear from the shallowness 
of the Ir edge (at about 2 keV) introduced by the mirror 
coatings. This edge feature is suppressed by the pileup of 
the high-count 1 keV photons, which correspondingly are 
depleted in the observed spectrum. 

No CIAO task exists to correct for pileup at present, 
and therefore fits to the core spectrum will be subject to 
a systematic error in the sense of underestimating the core 
flux density at low energies and underestimating the spectral 
index of the core, a (defined in the sense that the flux density 
S v cc v~ a% ). In addition, this early observation was pointed 
so that the core lies on the junction of two readout nodes of 
the ACIS-S3 chip. In such a location, it is difficult accurately 
to estimate the response matrix appropriate for the analysis 
of the spectrum, and a further systematic error is possible. 

Nevertheless, some indication of the properties of the 
core X-radiation can be deduced from the X-ray spectrum. 
A fair fit to the data is obtained by fitting to an absorbed 
power-law spectrum (Fig. ^ . The best-fitting spectrum has 
a spectral index a = 0.34 ±0.03 (corresponding to a photon 
index V — a + 1 = 1.34), consistent with the spectral index 
found by Garilli & Maccagni (1990) using EXOSAT, but 
flatter than the value a = 0.93 ± 0.03 reported by Hardcas- 
tle et al. (1999) based on the ROSATPSPC data. The effects 
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Figure 6. An expanded view of the core and jet region, with the 
same grey scale and contours as in Fig. [IJ The X-ray structure 
associated with the bright knot at the base of the radio jet is seen 
to lie roughly perpendicular to the radio axis. 



of pileup and inadequate calibration of the response matrix 
may cause a systematic error of about 0.1 on the measured 
power-law index. The corresponding best-fit 1 keV photon 
flux density, (1.03 ± 0.02) x 10~ 3 photons cm" 2 keV" 1 s" 1 
(0.68 ± 0.01 mjy), is likely to be a serious underestimate 
(and, indeed, is a factor « 3 below the flux density measured 
by the ROSAT PSPC). If account is taken of the Galactic 
column towards PKS 0521-365 (N H = 3.37 x 10 20 cm" 2 ; 
Elvis, Wilkes & Lockman 1989), no absorbing column asso- 
ciated with the AGN is required in this fit (Fig. ||): the 2a 
upper limit is A*h < 8 x 10 19 cm" 2 . Pileup and calibration 
uncertainties affect the shape of the core spectrum below the 
peak, but an intrinsic absorbing column significantly larger 
than this limit would strongly affect the cut-off energy in 
Fig. |i|, so that we believe that the limit on Ah quoted above 
is reliable, and indicates that we are observing the core emis- 
sion free from an absorbing column as large as found towards 
the cores of FRIRGs. 

No correction for background was needed in making 
these spectral fits: a nearby background region with the 
same area as the source extraction region contained only 
4 counts (about 3 of which arise from the wings of the re- 
sponse from the core), compared with the 7900 counts in the 
spectrum shown in Fig. ]d. 



3.2 The jet 

A close examination of the structure of the X-ray emission 
from the jet knot (knot A in the notation of Scarpa et al. 
1999) appears to show an elongated structure with its long 
axis lying in position angle (7 ±2)°. This is roughly perpen- 
dicular to the local ridge-line of the radio jet (see Fig. |6j), 
and differs significantly from the position angles of the chip 
pixelization. The elongation may therefore be real despite 
appearing exceptionally narrow, and potentially being sub- 
ject to dither effects and aliasing in the ACIS pixel assign- 
ments. Similar narrow, jet-crossing, features have been seen 
in other radio jets (e.g., 3C66B; Hardcastle, Birkinshaw & 
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Figure 7. The spectrum of a 2.2-arcsec diameter region cen- 
tred on the X-ray knot (upper panel) and the residuals from the 
fit (lower panel). The best fitting power-law model for the knot 
spectrum is shown as the solid line, and corresponds to energy in- 
dex a = 1.41 ± 0.57 (photon index T = 2.41), 1-keV flux density 
15±5 njy, and an absorbing column of (1.7±4.7) X 10 20 cm" 2 at 
the redshift of PKS 0521—365. The background from the wings of 
the point response of the core has not been subtracted in this spec- 
trum, and is simultaneously fitted. A good description of the back- 
ground is given by a power law with energy index a = 1.58 ±0.20 
and a 1-keV flux density about 45 per cent that of the jet knot. 
The quality of the fit is good: \ 2 = 13.9 with 23 degrees of free- 
dom in background and source combined. 



Worrall 2001b), and tend to lie just core- wards of bright jet 
radio knots. 

The emission from this X-ray knot is clearly contami- 
nated by counts from the wings of the response to the AGN. 
334 counts at energy < 10 keV are found from knot A, where 
the counts were extracted from a circular region with diam- 
eter 2.2 arcsec centred on the knot. A similar extraction 
region on the opposite side of the core was used for back- 
ground, since it was expected to be equally contaminated by 
counts from the wings of the core. This background region 
contains 92 counts. Neither the jet knot nor the background 
region show count rates high enough that pileup is a prob- 
lem. Thus about a quarter of the counts from the region 
centred on the jet knot can be attributed to contamination 
from the core, or emission from the thermal halo. After cor- 
recting for this background, we find that the knot is detected 
by Chandra at 0.031 ± 0.003 counts" 1 in 0.3 - 7.0 keV. 

The spectrum extracted from the jet region including 
background from the core and thermal halo is shown in Fig. m 
Sufficient counts exist for meaningful fits only in the en- 
ergy range 0.3 — 4.4 keV. The fits were performed using 
SHERPA's ability to fit source and background simultane- 
ously. The jet emission was modelled as a power-law, ab- 
sorbed by the Galactic column and a possible excess at the 
redshift of PKS 0521—365. The background was modelled as 
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Table 1. Jet knot flux density measurements 
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Figure 8. \ 2 contours for fits to the spectrum of the jet of 
PKS 0521— 365 as a function of the photon index of the fitted 
spectrum, rj Gt , and the intrinsic absorbing column in front of 
the jet, TVjj- Contours are drawn at 1, 2, and 3c uncertainty 
offsets from the best fit, which lies at \ 2 = 13.9, T = 2.41, 
A^h = 1.7 X 10 20 cm -2 . The best fit with no such absorption 
has x 2 = 14.0 (with one additional degree of freedom) and pho- 
ton index T = 2.25. 



a simple power-law (found to have a similar photon index 
to that of the jet knot; Fig. ^), which was expected to be 
a good representation of the superposition of the peripheral 
counts from the core of the the AGN and the thermal halo. 
The best fit (shown in Fig. |7|) involves an intrinsic ab- 



sorbing column of 1.7 x 10 cm" , but this is not significant 
(Fig. p|). If the additional column is set to zero, as suggested 
by the fit to the core spectrum, the jet energy spectral index 
changes from 1.41 ± 0.57 to 1.25 ± 0.30 (la errors with one 
interesting parameter). The significant change in the error 
on the spectral index arises from the strong correlation be- 
tween the fitted power-law index and the poorly-constrained 
absorbing column (Fig. ^). The value of \ 2 a t the best fit 
(13.9, with 23 degrees of freedom) is only slightly less than 
the value with no absorbing column (x 2 = 14.0, with 24 de- 
grees of freedom). None of the other fitted parameters vary 
by more than 10 per cent of their errors if it is assumed that 
there is no intrinsic absorbing column at PKS 0521—365. 

The 1-keV flux density of the jet knot, from the best fit 
without intrinsic absorption, is 14±3 njy, which corresponds 
to an 0.2 - 2.4 keV luminosity of 1.0 x 10 35 W if the emission 
is isotropic. The spectral energy distribution of this knot is 
shown in Fig. ^| where flux densities (given in Table are 
for a nominal 2.2-arcsec diameter aperture centred on the X- 
ray knot. The solid line drawn on Fig. ^ is calculated as the 
synchrotron emission from a population of electrons with a 
broken power-law energy distribution. 

On Fig. [] the radio spectrum appears steeper by about 
0.2 than the model, which has a radio spectral index of 0.58 



Frequency (Hz) 
4.9 x 10 9 
8.6 x 10 9 
1.5 x 10 10 
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60 ± 5 mjy 
45 ± 6 fiJy 
14 ± 3 njy 



Reference 
this work 
this work 
Scarpa et al. 
Scarpa et al. 
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(corresponding to an electron energy index p = 2.16), prob- 
ably because of a contribution to the lower- frequency flux 
densities from steeper-spectrum emission from the halo. The 
X-ray spectrum is significantly steeper, with a ~ 1.25. This 
extrapolates back to a good match to the R-band optical 
flux measured by Scarpa et al. (1999), and implies that the 
spectrum must break by Aa « 0.7 in the infra-red. For the 
equipartition magnetic field strength of 16 nT (estimated as- 
suming that the knot is moving non-relativistically, has an 
angular diameter of 0.5 arcsec, is completely filled by mag- 
netic fields and relativistic electrons of energies > 150 MeV, 
and that there is a negligible contribution to the total en- 
ergy from protons), this break lies at an electron energy of 
about 200 GeV. 

The amplitude Aa ~ 0.7 of the break in the spectrum 
corresponds to a change in electron energy power-law index 
Ap ~ 1.4. Such a large break is unexpected from simple 
synchrotron ageing arguments, and suggests that there is 
a change in efficiency of the electron acceleration process 
at about 200 GeV. In the shape of the inferred electron en- 
ergy distribution and the value of the equipartition magnetic 
field, this knot is very similar to that in 3C 66B. 

Also plotted on Fig. [] are the expected levels of inverse- 
Compton emission from the electrons scattering the mi- 
crowave background radiation and the jet's own synchrotron 
emission, assuming that the population of electrons in the 
jet is given by equipartition arguments. The turndown in 
the predicted inverse-Compton emission from scattering the 
microwave background radiation occurs because of the as- 
sumed lower energy limit of 150 MeV for the electron pop- 
ulation. Neither scattering process can be responsible for as 
much as 10 -3 of the X-ray emission seen under the assump- 
tion that the jet is non-relativistic. 

If the jet speed is highly relativistic, then the level of 
inverse-Compton emission from scattering the microwave 
background is higher than modelled in Fig. bl with a Lorentz 
factor ~ 30 being sufficient to produce the observed X-ray 
brightness. However, even in this case the predicted spec- 
trum would be too flat to match the steep X-ray spectrum 
measured by Chandra unless there is a substantial popu- 
lation of low-energy electrons (which must emit radio syn- 
chrotron radiation only at frequencies below 5 GHz). An ad- 
ditional photon source, for example beamed emission from 
the nucleus, can also add extra inverse-Compton X-ray emis- 
sion, but again the predicted spectrum would be a poor 
match to the data unless a low-energy electron population 
is present. 

Thus, since the X-ray spectrum produced by inverse- 
Compton X-ray emission would have a ps 0.6 (the same 
spectral index as the source photons, which are upscattered 
principally from the radio and mm- wave bands), while the 
observed spectrum has a = 1.25 ± 0.30, we conclude that 
the jet X-ray emission in PKS 0521— 365 arises from syn- 
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Figure 9. A synchrotron model roughly matching the radio, op- 
tical, and X-ray flux densities of the bright knot in the jet of 
PKS0521-365as given in Table [[ The R-band and 15-GHz data 
are taken from Scarpa et al. (1999), the other radio flux densi- 
ties are from our ATCA and archival VLA data, and the X-ray 
spectrum is shown with the 1-keV flux density and the permitted 
range of spectral indices assuming that the intrinsic absorbing col- 
umn is zero. The model requires a break in the infra-red of about 
Aa = 0.7, corresponding to an electron energy of about 200 GeV 
in the equipartition field. Neither the synchrotron self-Compton 
radiation (dash-dotted line) nor the inverse-Compton radiation 
produced by scattering the microwave background radiation (dot- 
ted line) makes a significant contribution to the observed X-ray 
flux of the jet. 

chrotron radiation, as it does in FRIRGs (Hardcastle et al. 
2001b; Worrall et al. 2001). The appearance of the spec- 
tral energy distribution from PKS 0521-365 (Fig. |§) also 
strongly resembles that seen in FRIRGs. 

3.3 The halo 

Hardcastle et al.'s (1999) analysis of the ROSAT HRI data 
implied a bright X-ray halo about PKS 0521— 365. At the 
best-fit halo brightness, and assuming that the temperature 
of gas responsible for the halo is « 1 keV we would have 
expected Chandra to detect about 10 times as many counts 
from the halo in Fig. ^ as were found in the ROSAT image. 
However, the radial profile (Fig. |^) shows a much fainter 
level of emission. At 10 arcsec from the core, where the ex- 
cess above the PRF wings is strongest, the halo contributes 
about 5 x 10 -5 count arcsec -2 s -1 . The earlier fits to the 
ROSAT HRI data showed a best-fit brightness at 10 arcsec 
of about 2 x 10 -5 HRI count arcsec -2 s - , which would sug- 
gest an ACIS brightness of 2 x 10 -4 count arcsec -2 s , four 
times brighter than the halo that appears in Fig. ^. 

However, the error on the brightness of the halo de- 
rived from the HRI image analysis was considerable, since 
excess counts above those expected for a bright point source 
were seen only in the angular range 10 — 30 arcsec, with the 
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Figure 10. The spectrum of the X-ray halo (upper panel) and the 
residuals from the fit (lower panel). The best fitting Raymond- 
Smith thermal model is shown as the solid line, and corresponds 
to a model with k-gT = 1.58 keV, abundance 0.06 solar, and no 
absorbing column above the Galactic level. Background emission 
is also present in this spectrum, and contributes about half the 
total counts (especially above 3 keV). This background is simul- 
taneously fitted, and can be well described by two power laws, 
with photon indices Y = 2.46 and —0.13. The quality of the fit is 
good (x 2 = 13.3 with 35 degrees of freedom in background and 
source) . 



wings dominating further out and the core response domi- 
nating close in. If a slightly altered HRI point response func- 
tion (based on more recent calibrations than were available 
at the time of the Hardcastle et al. analysis) is used, the 
best fit changes dramatically, from /3 — 0.90 to 0.66, core 
radius 8 to 3.6 arcsec, and halo brightness at 10 arcsec of 
2 x 10 -5 to 5 x 10 -5 HRI count arcsec -2 s -1 . While these 
results are even less consistent with the Chandra data, the 
strong change in the fit based on a small variation in the 
HRI calibration indicates that large systematic errors are 
possible in the HRI-derived halo parameters. Furthermore, 
when account is taken of the range of models allowed by 
the HRI data, the error on the 10-arcsec count rate is about 
50 per cent. 

If the radial profile in Fig. ^ is fitted using a standard 
point response for the core, and the central region is ex- 
cluded because of the misfit anticipated because of pileup, it 
is found that the best-fitting /3 model halo has structural pa- 
rameters = 0.90 and core radius 9.3 arcsec, strikingly sim- 
ilar to the values found by Hardcastle et al. (1999). Adopt- 
ing the new Chandra normalization of the halo, however, has 
the effect of reducing the halo brightness by about an order 
of magnitude from the value reported based on the ROSAT 
data: the implied total X-ray emission from the halo of about 
0.065 counts -1 in 0.3 — 7.0 keV, corresponds to a bolometric 
X-ray luminosity Lx ~ 3 x 10 35 W (using a temperature of 
1.5 keV for the halo gas), or 2 x 10 35 W in the 0.2 - 1.9 keV 
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Figure 11. \ 2 contours for fits to the spectrum of the halo of 
PKS 0521—365 as a function of the temperature of the fitted spec- 
trum, fcB^haloi an d the metal abundance of the gas, A^^. Con- 
tours are drawn at 1, 2, and 3<r uncertainty offsets from the best 
fit. 

band, about a factor 4 fainter than was measured by ROSAT 
(Hardcastle et al. 1999). 

A spectral analysis of the halo was attempted by ex- 
tracting the counts from a region between 5 and 20 arc- 
sec from the core, excluding the readout streaks, and back- 
ground from 20 to 27 arcsec from the core. The halo extrac- 
tion region contains 554 counts, while the background region 
contains 212 counts of which about 27 (based on the best- 
fitting structural model) arise from the outer halo. Roughly 
40 counts in the halo region arise from the wings of the re- 
sponse to the AGN. This complicated situation was dealt 
with by making a simultaneous fit to the halo and back- 
ground spectra, with the background represented by two 
power-law spectra: a flat component that models the parti- 
cle background, and a steep component that models the soft 
X-ray background. This representation of the background 
is consistent with the spectra seen in other regions of the 
S3 chip. A good fit (Fig. 0) is found with the halo emission 
arising from a (Raymond-Smith) thermal plasma with tem- 
perature feT = 1-6^0 4 keV, abundance 0.05ln'riR solar (1°" 
errors with 2 interesting parameters, see Fig. Ill]) , and 1-keV 
flux density 130 ± 20 njy. No additional absorbing column 
is needed. A thermal spectrum for the halo, with % 2 = 13.3 
and 35 degrees of freedom, is superior to a power-law fit 
(X 2 = 16.5 with 36 de grees of freedom) . Adding a power-law 
component to the halo model, to correct for counts from the 
wings of the AGN's PRF, has no significant effect on the 
goodness of fit or the best-fit parameters. 

The halo temperature and luminosity are similar to the 
values seen from the atmospheres of FRIRGs and lie close to 
the correlation between these quantities described by Wor- 
rall & Birkinshaw (2000). Thus the properties of the halo of 
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Figure 12. The hotspot of PKS 0521-365 from our 8.6-GHz 
ATCA image, with the same contours as Fig. |l| superimposed 
on the 0.3-7.0 keV Chandra image convolved to a resolution of 
0.75 arcsec. Note the close alignment of the radio hotspot with a 
relatively bright X-ray feature, which we believe to arise from 
inverse-Compton scattering of the synchrotron radiation by a 
near-equipartition plasma in the hotspot. 

PKS 0521— 365 are consistent with the unification hypoth- 
esis, that BL Lac sources and FRIRGs are from the same 
population seen at different orientations. 

3.4 The hotspot 

An excess of counts is seen within the 20 mjy contour of the 
radio map of Fig. |l|. This is better represented in the con- 
volved image of Fig. [l2| where a clear X-ray peak is associ- 
ated with the radio hotspot. We find 5 counts within 1 arcsec 
of the centre of the hotspot, while the local background im- 
plies that only 1 count would be expected by chance. This 
corresponds to a 3<r significant detection of the hotspot, al- 
though its flux, of about 5 x 10 -4 counts -1 in 0.3 — 7.0 keV, 
is poorly determined. 

The corresponding 1 keV flux density of the hotspot 
is about 0.4 njy. If we model it as an 0.4-arcsec diam- 
eter sphere (consistent with its extension on the 8.6-GHz 
ATCA map, Fig. |l]), and assume it to be completely filled 
with a uniform magnetic field and a relativistic electron- 
proton plasma whose energetics are dominated by the elec- 
trons, then assuming equipartition and electron energy lim- 
its 600 MeV - 400 GeV, the predicted synchrotron self- 
Compton X-ray flux density from the hotspot is about 
0.2 njy. 

While the significance of the detection is low, the rea- 
sonable agreement between the observed and predicted X- 
ray count rates from the hotspot is consistent with the 
hotspot being close to equipartition between its popula- 
tion of relativistic electrons and magnetic energy. Inhomo- 
geneities in the structure of the hotspot, the presence of 
electrons with different energy distributions, or small de- 
partures from equipartition are sufficient to match the ob- 
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served X-ray flux density more closely. Similar results have 
now been seen in a number of other radio hotspots (e.g., 
3C 123; Hardcastle, Birkinshaw & Worrall 2001a). 



4 CONCLUSIONS 

Our observation of PKS 0521— 365 has confirmed the exis- 
tence of an X-ray emitting halo around the AGN, but at a 
lower luminosity than suggested by the ROSAT HRI data. 
The brightness of the halo, as seen with Chandra is consis- 
tent with the halo brightnesses seen around FRIRGs, sup- 
porting the idea that this BL Lac is an FRIRG seen in a 
favourable orientation, and that BL Lac objects do not have 
unusual X-ray atmospheres. 

The detection of an X-ray bright region near the radio 
knot in the jet, though initially surprising, is again consistent 
with recent Chandra results on FRIRGs. Indeed, the spec- 
tral energy distribution of emission from the knot strongly 
resembles that from the inner jet region of 3C 66B (Hard- 
castle et al. 2001b) or M 87 (Bohringer et al. 2001), and 
can be fitted well as synchrotron emission with a broken 
power-law spectrum. However, the break Aa = 0.7 that we 
infer in the spectrum is larger than the value Aa = 0.5 
that would occur in a continuous-injection model for the 
relativistic electrons, and is also inconsistent with a single- 
injection aged spectrum, which would have an exponential 
cut-off above some maximum energy. The flattened appear- 
ance of the knot X-rays, their slight offset from the radio 
peak, and the analogy with FRIRGs, suggests that we are 
seeing evidence of energy-dependent electron acceleration. 

There is evidence that X-rays are being detected from 
the hotspot of PKS 0521— 365 at a level consistent with 
equipartition of relativistic electron and magnetic field en- 
ergy densities in this structure. This is consistent with recent 
results for the hotspots in many FR II radio galaxies. 

Further study of the intriguing emission structure in 
the radio jet demands better optical, near-IR, and sub-mm 
data to constrain the shape of the spectral energy distri- 
bution and confirm our inference about the location of the 
break, but could also benefit from a longer X-ray exposure 
to acquire better statistics on the shape of the emission re- 
gion in the X-ray. A longer X-ray exposure would also check 
our tentative detection of the hotspot. 

Clearly it would be of great interest to discover whether 
the results seen for PKS 0521— 365, which suggest that it is 
a beamed radio galaxy as required by unification schemes, 
are common to a sample of BL Lac objects rather than a 
peculiarity of this single source. PKS 0521— 365 is in some 
ways an unusual BL Lac. Its radio structure displays char- 
acteristics of both FR I and FR II sources oriented close to 
the line of sight. Although radio-selected (and hence a Ra- 
dio BL Lac), it has been described as a flat-spectrum radio 
quasar on the basis that it has a spectral energy distribution 
peaking at high energies (like the X-ray selected BL Lac ob- 
jects: more strictly, the high-energy peaked class rather than 
the low-energy peaked class) and a rather steep X-ray spec- 
trum (Donato et al. 2001). Furthermore, its spectrum has 
both broad and narrow emission lines above a strong blue 
continuum, plus emission from the stellar component of the 
underlying elliptical galaxy (Danziger et al. 1979), and so 
has some resemblance to a low-power, flat radio-spectrum, 



quasar, although these emission lines are relatively weak and 
the overall appearance of the spectrum is very similar to that 
of the well-known BL Lac 3C 371 (Angel & Stockman 1980), 
in which Chandra has also mapped jet X-ray emission (Pesce 
et al. 2001). It is therefore important that more low-redshift 
BL Lac objects selected to cover a range of properties, be 
well-imaged by Chandra to see whether similar X-ray jets 
and halos can be found. 
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